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Probabilistic determination of the body voltage in digital PD-SOI circuits and its 

application to static timing analysis 



Abstract 

We describe a technique for using knowledge of the average 
switching behavior of partially-depleted silicon-on-insulator 
(PD-SOI) circuits to reduce delay uncertainty produced by 
the floating-body effect. The approach is based on a unique 
state diagram abstraction of the PD-SOI FET which cap- 
tures all of the essential device physics, A simple analytic 
model of the body voltage is derived and calibrated based 
on Monte Carlo simulations. This model is used within the 
context of a prototype transistor-level static timing anal- 
ysis engine. Results are presented which demonstrate the 
dramatic reduction in delay uncertainty possible with this 
technique. 

r 

1 Introduction 

Silicon-on-insulator (SOI) technology has long found niche 
applications for radiation- hardened or high- voltage integrated 
circuits. Recently, however, SOI has emerged as a technol- 
ogy for high-performance, low-power deep-submicron digital 
integrated circuits[l]. For digital applications, fully-depleted 
devices have been largely abandoned in favor of partially- 
depleted technology, because of the difficulty in control- 
ling the threshold voltage of fully-depleted thin-film tran- 
sistors. Partially-depleted SOI (PD-SOI) has two main ad- 
vantages for digital applications: the reduction of parasitic 
source-drain depletion capacitances and the reduction of the 
reverse-body effect in stack structures and pass- transistor 
logic. 

At the device and circuit level, however, the floating 
body effect in partially-depleted SOI (PD-SOI) poses major 
challenges in the successful use of this technology. There can 
be large "uncertainties" in the body potential based on past 
switching activities. For many circuits, the design margin- 
ing required to protect against this uncertainty erodes all 
of the potential performance advantage under nominal op- 
eration. In addition, for many circuit styles in which noise 
margin is strongly determined by threshold voltage (e. g. 
dynamic circuits), considerable overdesign for noise can also 
result from conservative body-voltage margining. 

Previous work on PD-SOI has focussed on device issues[2] 
or delay effects due to the floating-body effect evident for 
particular circuits under periodic stimulus[3, 4] (pulse stretch- 
ing, frequency-dependent delay time). No technique has yet 
been developed for verifying floating-body effect on timing 
systematically across circuit styles. In this paper, we con- 
sider characterizing the body voltage uncertainty of PD-SOI 
devices using probabilistic knowledge of switching activity 
and applying this to transistor-level static timing analysis. 
In our approach, we analyze each FET of each circuit, de- 
termining the minimum (min) and maximum (max) body 
voltage under both "initial-condition" and "steady-state" 
operation. By initial-condition operation, we mean that the 



circuits are assumed to be quiescent for a long time and 
then undergo one cycle of switching activity to produce the 
min/max body voltage. Steady-state body voltage variation 
takes into account the average switching behavior of the cir- 
cuit and is always considerably less conservative than the 
initial condition values. This technique formalizes the ob- 
servation that one can expect very different behavior from 
a functional block on-chip that is quiescent for a long pe- 
riod of time as compared to one that has been under steady 
switching activity. 

In this paper, we work with BSIM3SOI[5] model pa- 
rameters for an IBM partially-depleted SOI technology de- 
scribed elsewhere[6]. Devices have a 0.25/ira effective chan- 
nel length, 5 — nm gate oxide, 350 — nm back oxide, and 
140 — nm thin silicon film. The supply voltage is 2.5V. While 
the detailed results we present here apply to this technol- 
ogy, the techniques are generally applicable to any PD-SOI 
technology. 

In Section 2, we describe a state-diagram abstraction of 
the PD-SOI FET which simplifies the device physics deter- 
mining the body voltage. In Section 3, we describe how we 
analyze this state diagram to determine the initial-condition 
and steady-state min/max body voltages. Section 4 de- 
scribes a prototype transistor-level static timing analysis 
engine which incorporates these body voltage characteriza- 
tions. Some results are presented in Section 5. Section 6 
concludes and offers direction for future work. 

2 State diagram view of body interactions 

The body potential of a PD-SOI FET is determined by ca- 
pacitive coupling of the body to the gate, source, and drain 
and by diode currents at the source-body and drain-body 
junctions. To model the switching history determining the 
current body voltage of a particular device, we use the state 
diagram abstraction shown in Figure 1 for an nFET. Ar- 
rows indicate possible state transitions produced by switch- 
ing event, while the states themselves represent valid logic 
conditions for an FET at the end of a cycle or during haz- 
ards that occur transiently within a cycle. The states 1, 2, 
5a, and 5b are distinguished in Figure 1 because transitions 
between these states propagate timing delays. States 5a and 
5b can usually be treated equivalently as state 5. The state 
diagram of the pFET is the "dual" of Figure 1, in which the 
gate is high rather than low in states 3, 4, and 5; and low 
rather than high in states 1 and 2. 

In analyzing the device physics determining the body 
potential, it is convenient to distinguish "slow" and "fast" 
processes. Fast processes can change the body potential on 
time scales on the order of or less than the cycle time, while 
slow processes require time scales much longer than the cycle 
time to affect the body voltage. There are two fast mech- 
anisms at work: switching transitions on the gate, source, 
or drain which are capacitively coupled to the body, and 
forward-bias diode currents across source-body and drain- 
body functions with voltages exceeding the diode turn-on 




Figure 1: State diagram for a PD-SOI nFET. 

voltage. The latter process results in a fast discharge of the 
body through the forward-biased junction. The coupling 
"kicks" associated with the first mechanism occur for each 
transition in Figure 1 and are completely reversible on "fast" 
time scales; that is, if one begins in state 1 and traverses the 
state diagram, returning to state 1 on a time scale compa- 
rable or less than the cycle time, the body voltage on return 
will be the same as the body voltage. The reversibility of 
the kicks is, of course, dependent on not triggering the irre- 
versible discharge mechanism. In practice, the fast discharge 
mechanism is only triggered on activation of a circuit after 
a long period of quiescence. Once the discharge happens, it 
is not triggered again as long as the circuit is under steady 
switching activity. 

The slow processes involve charging or discharging the 
body through reverse-biased or very weakly forward- biased 
diode junctions. These leakage currents give each state a 
slow relaxation to a target dc body voltage, which is indi- 
cated for the nFET case in Figure 1 for our target tech- 
nology. We note that the dc body voltage of state 5 is 
determined by the steady-state balance between a weakly 
forward-biased junction and a reverse-biased junction. For 
the pFET case, the dc voltage of state 5 in our target tech- 
nology is 1.73 V. 

3 Determining the body potential 

Important to determining the steady-state body voltage is 
characterizing the average switching behavior of the circuits. 
To do this, we appeal to stochastic techniques to deter- 
mine the average fraction of time an FET spends in each of 
states 1 through 5. This is determined by two factors: the 
fraction of cycles an FET ends the cycle in a given state; 
and the amount of time the FET spends in each state with 
each cycle-to-cycle transition (including the possibility of 
hazards). We also find, based on simulation results, that 
to maximize (minimize) the body potential, we must favor 
those states with high (low) source-drain over those with 
low (high). In addition, although less importantly, we favor 
state 4(1) over state 2(3) for maximizing (minimizing) the 
body voltage for the nFET. For the pFET, we favor state 
1(4) over state 3(2) for maximizing (minimizing) the body 
voltage. We expect these observations to extend generally 
to other PD-SOI technologies. 

Each FET in our analysis is characterized by a set of 
arrival times and signal probabilities. For the nFET, these 
are: J 

• P{G) = probability that at the end of a cycle the gate 
is high 



• P(D\G) = probability that at the end of a cycle the 
drain is driven high given that the gate is low 

• P(D\G) = probability that at the end of a cycle the 
drain is driven low given that the gate is low 

• P(S\G) = P(D\G) = probability that at the end of a 
cycle the drain (and source) are driven high given that 
the gate is high 

• P(S\G) = P(D\G) = probability that at the end of a 
cycle that the drain (and source) are driven low given 
that the gate is low 
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times (rising and falling) for the source, gate, and drain 

By driven low (high), we mean that there is a path to ground 
(supply). The details of how these probabilities and arrival 
times_are determined * s described in Section 4. We note that 
P(D\G) + P(D\G) < 1, for example, but does not have to 
sum to 1, since the conditional signal probabilities do not 
cover the case that the drain is floating. 

Given these probabilities, the state diagram in Figure 
1 represents a Markov process with two six-by-six transi- 
tion matrices A (with one transition per machine cycle) 
corresponding to the min and max cases: {P m ax)k+i = 
Amax{PLax)k and (P min )k+i = A min (P min ) k , Each col- 
umn of A ma;r/min is given by 

/ max / min max f min max /mm max I min max /min max I min\T 
(Pi P2 V Z Pi P5a Ptb ) • 

Pi is the probability of making a transition to state i. Pi 
is the probability of being in state i. For example, p max = 
(1 - P(G))P(S\G)P(D\G) while pf n = (l - P{G)){1 - 
P(S\G))(1 - P(D\G)). The max (min) case assumes that 
the floating node condition goes high (low). Diagonaliz- 
ing A (trivially) and finding the eigenvector associated with 
eigenvalue 1 (normalized so that the sum of the elements 
of the vector is 1) gives the steady-state values of the Pi. 
From these probabilities, one can calculate a set of thirty- 
six transition probabilities Pi—> that define the probability 
that a FET is current in state i transitioning to state j: 

For each transition, we determine the fraction of the cycle 
time (tl~* k ) that can be spent in each state i as part of 
the transition j : — j- k to maximize (or minimize) the body 
voltage. For example, if the nFET is transitioning from state 
3 to state 3, then to maximize the body voltage: 

t>- 3 = Max(D'tf h SftZ) - Mtn(D%:?,SZ;») 

j3 — * 3 i j 3 — <■ 3 

l 3 — f'cycle M 

where tcycie is the cycle time and the source and drain are 
not constant 0 and, therefore, can execute such a hazard 
during a cycle. From these times, we can calculate an effec- 
tive amount of time (t***) on the average per cycle that the 
FET is in the state i: 

We now seek an analytic expression that allows us to deter- 
mine the steady-state voltage from these t^* . 

If we imagine an average cycle broken into slices by t*** , 
the body voltage in this average cycle will appear as in Fig- 
ure 2. The discontinuities between the slices occur as a result 
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Figure 2: An average cycle. 

of the capacitive coupling kicks and are characterized by the 
voltage differences di from the reference level (in this case, 
for the nFET, the body voltage of state 1). In addition, 
each state is relaxing (very slowly, which is why any voltage 
change during the time slice is imperceptible in Figure 2) 
to the target dc values noted in Figure 1 and here denoted 
as Si. The time constants for this relaxation are denoted as 
n. From this simple picture, one can relate the body volt- 
age at the end of the cycle v n +i to the body voltage at the 
beginning of the cycle v n by: 

«»+! = v n e£*=r ulTi + £( 4i - «,■)(! - e-'i'^e^U '' /ry 



steady-state. The steady-state value is independent of the 
the initial condition. 
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Figure 3: Body voltage transient. 



The steady-state solution of this difference equation (in the 
approximatioin that the r» are much greater than the t*^) 
is given by: 

E,(Wn)(^-rf.) 

We treat the relaxation time constants n as fitting parame- 
ters that we determine to match the steady-state body volt- 
ages that result from Monte Carlo simulations of the state 
diagram in Figure 1. One such simulation is shown in Fig- 
ure 3 for state 1 of the nFET and state 2 of the pFET. In 
this case the cycle time is 5 ns, all of the early arrival times 
are set to 1 ns, all of the late arrival times are set to 4 ns, 
and all of the signal probabilities are 0.5. In Figure 4, to 
demonstrate the efficacy of the analytic body voltage model, 
we compare it with Monte Carlo simulation results for an 
nFET (min and max). The conditions are the same as for 
Figure 3, except the signal probabilility of the gate is varied 
from 0 to 1. 

In the case that the circuits can be quiescent for a long 
period of time before undergoing switching activity, one 
must calculate the minimum and maximum body initial- 
condition body voltages for states 1, 2, and 5. This calcu- 
lation is done assuming that the device achieves a dc bias 
at one state and a capacitive kick (on a transition) into 
the target state. Based on simulation results for our target 
technology, we find that to maximize the body voltage, the 
initial state should be chosen with the priority: 4, 2, 5, 3, 
1. Lower priority initial states are chosen only if high pri- 
ority ones are probabilistically inaccessible. For the pFET, 
the initial state priority for the max case is: 2, 4, 5, 1, 3. 
The min case priority is the inverse of the max case priority. 
These initial conditions are noted for the example of Figure 
3. The rapid progression from this state near the start of 
the simulation is due to forward-bias discharge of the body 
(a fast process). Slow processes alone then determine the 




Figure 4: Body voltage fit. 



4 Static timing analysis 

We have incorporated this body voltage characterization 
into a prototype static transistor-level timing analysis engine 
utilizing a breadth-first-search algorithm. The design is par- 
titioned into channel-connected components (CCC) for anal- 
ysis, as is traditionally done in static transistor-level tools. 
The primary distinguishing feature of this timing analysis 
engine is that it propagates switching probabilities using as- 
sumptions of spatial and temporal independence, borrowing 
from similar techniques in static power analysis[7]. Once 
the signal probabilities and times are known at the inputs 



of a CCC, these probabilities are translated into FET sig- 
nal probability and arrival time properties for body-voltage 
calculation using an exhaustive path trace through the tran- 
sistors of the CCC from the target node to supply or ground. 

5 Results 

Because of the length constraints of this paper, we present 
only two timing examples. The first is a 15-stage inverter 
chain. In Figure 5(a), we present the steady-state delay 
(for both rising and falling transitions) for each stage in 
the case that the switching probability at the input is 0.5. 
The minimum and maximum delays are virtually identical. 
We constrast this with Figure 5(b), in which we present the 
steady-state delay for the case that the switching probability 
at the input is 0.9. In this case, a wide disparity in the rising 
and falling delays opens up. In the case of the inverter, the 
only states accessible to the nFET are 1 and 5a f and for the 
pFET, 2 and 5b. For the inverters with input probabilities of 
0.9, the steady-state body voltages of the pFET and nFET 
are pulled below their 0.5 values. This makes the nFET 
faster and pFET slower. The opposite applied to the invert- 
ers with input probabilies of 0.1 in the chain. The result is 
the well-known "pulse-stretching" phenomenon [4], now ac- 
cessible within the static timing analysis environment. 

We next consider an static eight- bit ripple carry adder 
using the full-adder implementation of Reference [8]. We 
consider the critical path of this circuit, the carry chain. 
In Figure 6(a), we show the stage delay using the initial- 
condition values of the body voltage. The minimum and 
maximum rising and falling delays are noted. The very 
large delay uncertainty is dramatically reduced by introduc- 
ing knowledge that the circuit is under steady switching ac- 
tivity, as shown in Figure 6(b), which uses the steady-state 
analysis and assumes that all of the inputs have signal prob- 
abilities of 0.5. There is very little difference between the 
min and max cases in steady-state (very little delay uncer- 
tainty). 
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Figure 5: Inverter chain. 




Figure 6: Carry chain. 



6 Conclusions 
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We have demonstrated an efficient technique for bounding 
delay uncertainty in digital PD-SOI circuits using knowledge 
of the average switching behavior of the circuits, incorporat- 
ing "hysteric" effects into the context of static timing anal- 
ysis. Designers will naturally be very uncomfortable at first 
with applying stochastic techniques to something as "life- 
and-death" as timing or noise analysis. What they will find, 
however, is that the body- voltage uncertainty can be signif- 
icantly reduced with fairly conservative assumptions about 
switching behavior. In many cases, a normally inactive 
block may be periodically stimulated to keep it "primed" 
so that when it is eventually exercised, it has more tightly 
controlled body voltage variation. Determining the neces- 
sary frequency and nature of this pattern and how it could 
be incorporated practically into the design remains an open 
research question. 
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